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LASER PROPAGATION CODE STUDY 



by 

Edward B. Rockower 
I. INTRODUCTION and SUMMARY 

During the course o-f this study a number of laser pr opaq^^ L i or< 
codes have been assessed as to their suitability tor fnodeiinq 
Army High Energy Laser weapons used in an anti— sensor mode. 
Because the Army battle-field scenario requires non 1 i ne<.^.r 
beam propagation throi.iqh turbulent atmospheric conditions. 

1 1 1 c i Li d 1 n q s (x\ a f : e a n d d Li s t , t li e s e aspects ha v Oz? r e q li i r e d h i i .j 
received) tkie greatest modeling effort (see, e.q. LiJ-. 

However, other important feattires of High Energy Laser' weapoii. 
are the rectangular aperttire and large (possiblv Lip to vu'.; 
central obscuration necessitated b/ the unstable resonator 
designs ccHumon to most hiqh energy lasers- 

The adequacy oi the codes in modeling these, and other, 
■rspects or oper a I j. oru-t 1 lasers has required a I. i dat i or i . r e 
report will identify some apparent problems with current 
rnocielinq in the codes, propose an inteiiin f 3 •■ " , wliere i'Cs-: 1 

an‘j outline an approach to r ecoinrnended ^Lirthcr oor r , on 

our stLidy of the i ssLies identified. 

Two msi .ior cateqnr'ies of code ha /e been i I'rves 1 1 Lia to* i m ou> 
study : 
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1. Fundamental wave optics codes; these codes start 
from the basic microscopic laws of E-M radiation and implement 
various phenomenol ogi cal models of the atmospheric turbulence 
and particulates. A number of schemes for numerically 
integrating the resulting approK i mat i ons to the coupled partial 
differential equations have been mechanized on computer t.see 
Table I for a listing of these and other codes). 

2- Moderate accuracy system-level scaling law 
codes; these codes start from phenomenol oqi catl and anailyticai 
approM i mat 1 ons to the E-M wave equations, and attempt to 
represent the details of propagation through the atmosphere wi l.n 
a few parameters (e.q. a phase integral; based on integrated 
properties of the atmospheric conditions, light intensity, wa e 
length, etc. These codes are based on experimentally and 
ruamerically (from the fundamental wave optics codes; derived 
data bases, 

fi large portion of the analysis reported here is rela.ti v'e to 
two baseline codes, one from each at the above categories. 
the first category is the 4-D code developed by Joe Hleck and 
J i rn Morris at Lawrence Livermore IMational Labor atory i!2 J . in the 

s£?cond, scaling law, category is the bKLPRG code developed Dy 
Harold Breaux of the Ballistic Research Laboratory, a version or 
which IS contained in the HELAWS code which models ermy high 
energy anti -sensor laser weaponsC3J. During the course of this 



study, this code has been modi-fied slightly while installing it 
on the Naval Postgraduate School (NFS) IBM 3033 computer svstem 
in order to make it compatible with the CHS operatinq system and 
IBM Fortran 77. 

A relatively large number o-f computer codes have been 
developed at various institutions in order to model high enerq 
laser propagation. Unfortunately, the capabilities and 
limitations at these codes are not alwavs immediatel / appa.reni- 
from a perusail of their documentation. Some o+ the 1 i t a c i orii 
are common to all of the codes; e.q., the code may onl / hanoie 
circulai' laser apertures. This limitatiori has. noi stoppeu c;i^ 
application the codes outside their range of- /ali«Jit^- 

limitations are those, such as assuming a vacuum iri whjcri ’ihr- 
laser propagates, which clearly preclude use of sucJ'. (_odes orher 
than in a Space Warfare scenario. In order to attempt to male a 
pr e 1 imi nar V assessment of what s available we have relied 
heavj.lv on teleplione con ver sat i ons with engineers and anaj £t = 
a t. V a r i o u. s i I'l s 1 1 1 u L i o n s r ta g a r d i n g i: heir’ c o « i e s . i r. i d d i • . i • i , 
Significant source of information has been a survey :>t Cc- 
ex I s 1 1 nq wave’-op tics, seal i nq 1 aw , and s l mp L i ^ i ad qeometr i 
Laser’ pr opaqai. i on codes vsee Fable 1; carried out b\ L’r . •Mune- 
P. Reilly and co-workers at W. J. bchafer Associaiies, Inc. in 



1979 L41. 



TABLE 1: PROPAGATION CODES 



CODE 


TYPE* 


DEVELOPED BY 


^ 

EXERCISED BY 


A PM 


SC 


HEL Systems Project Office 


Army 


BREAUX 


SC 


Ballistics Research Laboratory 


Army i 


COMBO 


SC 


Air Force Weapons Laboratory 


Air Force 


EAPM 


SC 


Charles Stark Draper Laboraory 


Charles Stark Draper Lab. 


ESP III 


SC 


United Technologies Research Center 


Air Force 


ESP IIIA 


SC 


United Technologies Research Center 


Air Force/Army 


ESP IV 


SC 


United Technologies Research Center 


Air Force 


GEBHARDT 


SC 


Science Applications, Inc. 


Army 


GUTS 


SC 


Air Force Weapons Laboratory 


Air Force 


HELP (PROPMD) 


WO 


Air Force Weapons Laboratory 


Air Force ^ 




wo 


Far Field, Incorporated 


Far Field, Incorporated 

_ . . 1 


LASE 


SL 


Science Applications, Inc. 


Navy ! 


IJ^SNEX 


GO 


Lawrence Livermore Laboratory 


Lawrence Livermore Lab. j 


LL/SL 


SL 


Lincoln Laboratory 


Lincoln Laboratory j 


LL/WO 


WO 


Lincoln Laboratory 


■ ' ■ 1 , -J 

Lincoln Laboratory | 


LLL/4D 


wo 


Lawrence Livermore Laboratory 


Lawrence Livermore Lab. 

1 


MPLAW 


SL 


Naval Research Laboratory 


Navy [ 


NOLEC 


SC 


Naval Ordnance Laboratory 


f 

Navy ! 


NRL/CHM 


wo 


Naval Research Lab/SAI 


Navy I 


NRL/JL 


sc 


Naval Research Laboratory 


Navy j 


PHILLIPS/SL 


SL 


Science Applications, Inc. 


Navy 


PSM 


WO 


Charles Stark Draper Laboratory 


csDL : 


SAICOM 


SC 


Science Applications, Inc. 


Navy 


SSPARAMA 


WO 


Naval Research Laboratory 


Navy 


UTRC/WO 


wo 


United Technologies Research Center 


United Technologies Resear 
Center 


ZAPM 


SC 


W. J. Schafer Associates, Inc. 


W. J, Schafer Associates, 



* WO - Wave Optics Code SC - Simplified Code 

SL - Scaling Law GO - Geometric Optics/Hydrodynamic Code 

I 
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II. APERTURE SHAPE PROBLEMS 



In con ver sat i ons with analysts and programmers regarding 
their codes and in study ot the documention on other codes, 
where available, it is apparent that virtually all non-wave- 
optics codes -force the user to model his laser as having a 
circul^xr beam profile at the laser aperture. As mentioned 
above, most high energy chemical lasers have rectangular 
apertures- In addition, there is often a central or non-CET‘rr. ^ 
obscuration of the beam profile at the laser aperture, caus-io o 
on e rn i r r or of the un s t able r eson a t or c on f i g u r a 1 1 on - i n 
attempting to exercise such a propa^qation code to si-nui-u:'? 
rectangular aperture lasers, it is reasonable to tr> to fna : r l i i 
both the same output power as well as output beam inr.ensit- s 
for the real laser. If the former condition is met, then the 
latter is equivalent to requiring that the area oi the c i r'cu 1 -r 
aperture (with or without a central obscur-at i on ^ be the sa.r.e 
f o r t hi e r o a J. I a e r - 

Thie following analysis is an attempt to estimate the 
seriousness of the limitation to ci r'cu I ar apertures and 1 c 3.0 
t o u.< r o p u sal f o i- a n j. 1 1 1 e r i m f i x . F i r i a 1 l '/ , jv e e u l:i g e s t 
passible course of further work on this problem. 

Pirt asv inptot 1 c approx i mat i on f or mu.l a f or ost i ma 1 1 r*.g chc- 
fraction of encircled energy within a given radiLis h,as been 
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derived for imaging systems with oddly shaped aper tures C 5 3 . 
Modulation Transfer Function (MTF) techniques were applied to a 
uniformly illuminated aperture having arbitrary shape and 
obscurations, resulting in the following formula, valid for 
asymptotically large values of r, 

E(r) = l-AfR/(2jrV) 



where 



A 

f 

r 

K 

E 



= laser wavelength 
= effective focal length 
= radial dimension in the focal plane 
= HEL aperture per i met er-to— area ration 
= normalized encircled energy, i.e. the 



fraction of the energy transmitted by the aperture that falls 
within a circle of radius r about the geometrical focal point. 

The main feature of this result for our purposes is that, 
apart tram the laser wavelength and system focal distance, the 
fraction of encircled energy depends only on the ratio R/r, 
where R is the ratio (per i meter /area) for the imaging system 
aperture, and r is the radial dimension in the focal plane 
within which one wishes to determine the fraction of total 
energy. Asymptotically, the fraction of encircled energy is 
independent of details of the shape, apart from the value of R. 
It IS assumed that the laser beam propagates in a linear medium, 
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with no etfects ot atmospheric f 1 ucfew-iat i ons , wind, etc. and that 
the laser aperture is uniformly illuminated. 

Analyzing this result, we conclude that, at least for the 
restricted conditions under which it was derived, the radius Kr i 
encircling a given fraction of energy increases linear! v with R. 
In other words, a beam from an aperture with double the -alue uf 
R will, asymptotically, spread twice as much, from di r tract ion . 
Hence, again subject to the limitations of app 1 i cab i 1 1 1 / , the 
Hipread of a beam ir» free space propagation will be pr opor 1 1 .,.n - 1 
to the ratio of the lasers perirneter to its area. 

r It e r e i s a o rri e w hi a I: s 1 1 ri liar e f i t- c L t h h i. s a l i r \ tu . i 

ccjnsequence of imperfections of laser wavefront gjj ^ I i i. a*- tne 

point where." the beam leaves the 1 aeer aperture- l-iha t o -fr-r •. ru? 
cause of the degraded beam quality ‘phase front nistor tioits* , 
whE'ther f rorn i nhornoqenei t i es i n the 1 asi nq med i urn , rni r r or or 
lens impiarf ections , etc., the effects on beam propauafion are 
rep)resented with a parameter known as the “beam qucxlity" • hh* . 

Till? value of M ls a Iwa/s qr easier ^ tliaii or' equal i:o I aiKi i _ .i n 

as “times dirfr action linn ted"; i - , the beam sprearj is il r i — 

the di f f r'act 1 on limited rate of beam di verqonce •: pr'opor 1 1 ono i tj 
la(iibda/D>- Comparing the two res'..'. Lts for beam spr'Cac , vjc 
identify a possible method of compensat i nq , at least pai'ciall., 
for rectangular apertures E^nd var'ious sizes and types e . a . non- 
central) of obscuration. Our conclusion is, for ^ree sp^ce 
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propaqation, that the asymptotic spread of lasers emanatinq from 

% 

two apertures having the same value of the product li*R will be 
the same. Hence, a better approximation to a rectangular 
aperture, possibly with a central obscuration, than just findinq 
the circular aperture with the same total area, is to also 
change the vailue of beam quality, M, so that the relation 

HR = H'R' 



is preserved- 

For example, we model a circular 

central obscuration with another circular 
area, output power, and with beam quality 
These two requirements lead to. 

Cl r cul ar f^per tur es ; 



aper tur e 
aper tur e 
given bv 



wi th a 
of the same 
M ' = H (R/R ) . 



M' = M n + /F)//(l - F) 
D' = D /(I - F) 

where F is the fraction obscuration. 



We have also derived similar equations for rectangular apertures 
with and witliout obscur at i ons . The results for the latter are 

presented in graphical form in Figure 1. Each curve in that 
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BEAM QUALITY QUALITY 



BEAM QUALITY VS RATIO OF RECTANGLE 



Figure 




RATIO OF APERTURE DIMENSIONS 



BEAM quality VS FRACTION CSSCURATiON 




1 . 



Beam qual i ty 
and Fraction 



as a -function o-f 
obscuration (F) . 



Length/Width ratio 



(L/W) 
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figure gives the beam quality of an ''equivalent” uniform 
circular aperture with the same area as the rectangular aperture 
whose length to width ratio (L/W) is shown to the right of the 
curve. The curve for a circular aperture with obscuration F is 
also shown - 

It remains to determine whether the above results have any 
relevance or usefulness for nonlinear laser propaqation within a 
turbul ent atmosphere. 

In order to estimate the magnitude of the problem and to 
determine whether the similarity tr ansf or mat i on derived above 
ifiay be useful we have exercised the HELAWS code t cor 1 1 a. i n i ng 
ERLPRG; with the following specific questions in mind: 

1) How adequate is the remedy of simply rep resen t l iiq 
rectangular aperture by using a circular one of equal area-* 

2) How adequate is it to represent a centrally obscured 
a p e r “ t u r~ e by ( n e a n s of a c ire u 1 a r " a p e r t u re irvi t h t In e s a m e a v r “ a 1 i 
area, power, and (therefore) intensity'? 

3) Gi'/en that we would like to try to represent dif;erent 
aperture shapes and obscurations with our similarity 

tr ansf or mat 1 on , how good a fix can be obtained in the presence 
^'f 1 i near atmospheric effects such as turbulence and wi nd ' 
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4) How good is our tr ansf or mat i on in face of nonlinear 



e-ffects, i.e. thermal blooming? 

5) How can the existing codes be improved, possibl with 
an extension of the similarity tr ansf or mat i on , or otherwise. 

With reference to the above questions, we now loot at the 
following 3 figures- The data were generated with the output 
from HELAWS, using variations of the base case parameters si-fc:a->n 
in Table 2- 

On each of the figures, curve \1) represents the resulLi:: .-.h' 

reference uniformly illuminated circular aperture of di A(neter- • 
one meter, bearri quality equal to one. Curve *2) preset. 
results for a uniformly illuminated circular aperture with 3. LU 
central obscuration but with the same area and beam qiu.ality 
(i.e. M - 1; as the base case. Since F - .i, one or olu" 
pr evi ou.s t or mul ae y i el ds 0 - 1 . - Fi nal 1 y , cur e • > 

presold: 5 thie results for a unifornH. > illuminated circular 
aperture mo ob acur a<t i on ) , equal area as before, buL with bcani 
quality iialculated from our formula to qive the same -^rec* 3‘j_-«ce 
cisyriipt ot 1 c: hcr-mi spread as the aperture with ItJ"'. ob scur st i ._-n ■ i . c 
curve Z) c Our formula presented above yields M -•= i . 3c . ihis 
can also, £ipprox i rnatel y , be read off of figure 4 rmiicri 3i»-jws .I 



TABLE 2: HELAWS CODE INPUT PARAMETERS 



LASER PARAMETERS 

LASER CODE R1-B3 

OPERATING MODE RP 

WAVELENGTH 10.6 MICRONS 

BEAM QUALITY (TIMES DIFFRACTION LIMITED) 1.39 X DL 

APERTURE MIRROR DIAMETER 1.000 M 

TOTAL SYS. JITTER (1 SIGMA) 10.0 URAD 

VARIABLE FOCAL RANGE EQUAL TO TGT RANGE 

AIMPOINT OFFSET: X COORDINATE 0.00 M 

AIMPOINT OFFSET; Y COORDINATE 0.00 M 

STD. DEV. OF AIM BIAS ERROR (1 SIGMA) • . 5.0 URAD 

ENERGY PER PULSE 10.0 KJ 

PULSE REPETITION RATE 5.0 P/SEC 

PULSE DURATION (OR WIDTH) 20.0 USEC 

ATMOSPHERIC PARAMETERS 

SCENARIO LOCATION EUROPE 

SCENARIO LONGITUDE 9.0 DEG 

SCENARIO LATITUDE 50.0 DEG 

TIME OF YEAR OCT 

TIME OF DAY 10.00 HRS 

VISIBILITY RANGE 7.00 KM 

AMBIENT TEMPERATURE 10.0 DEG C 

RELATIVE HUMIDITY 85.0 

ATMOSPHERIC PRESSURE 985.0 MB 

TURBULENCE LEVEL MODERATE 

REFRACTIVE INDEX STRUCTURE CONST (AT IM) 8.40E-14 M**(-2/3) 

SCALING OPTIONS EMPLOYED YES 

MAGNITUDE OF WIND VELOCITY AT REF HEIGHT 1.0 M/S 

REFERENCE HEIGHT FOR WIND VEL ALTITUDE SCALING 1.0 M 

WIND DIRECTION ANGLE (TO S0UTH=0 DEG, TO WEST=90 DEG) 220. M 

ENGAGEMENT PARAMETERS 

INITIAL TARGET RANGE 4000. M 

MAXIMUM TARGET RANGE 4500. M 

RANGE INCREMENT 100. M 

ALTITUDE OF LASER APERTURE 2.8 M 

ALTITUDE OF TARGET 1.5 M 
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LASER PROPAGATION CODE STUDY 



Figure 3. Unbloomed peak fluence vs. range 
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LASER PROPAGATION CODE STUDY 



required beam quality o-f an equivalent circular aperture laser 

for various shapes and values of F. 

First, let's look at figures 2 and 3, showing the unb loomed 
1/e spot diameter and unbloomed 1/e peak fluence, respectively. 
Comparing curves (1) and (2) shows that simply trying to 
represent a 10-/. obscurat i on ' wi th a circular aperture of equal 
area would lead to gross errors. We now appeal to the encircled 
theorem and note that it is easy to show that there is a 
rectangular aperture (L/W approK. = 1.9) with the same area and 

value of R as the 10V. obscured circular aperture. Hence, bv the 
encircled energy theorem the rectangle will have the same 
asymptotic beam spread as will the obscured circle. Clearlv. if 
one tries to represent such a rectangle with the unobst-ui ed 
circular aperture (curve 1), then the same gross errors will 
result, even without thermal blooming (there's certainly nu 
reason to believe that the nonlinear thermal biooming will help 

matters here) . 

Still looking at figures 2 and 3, now note curve > -■; whi._.i 
corresponds to an unobscured circular aperture with che sume 
area as the other cases but with beam quality M set by the 
requirement that H*R matches that of the 107. obscuration case. 
We see that the match is eiscellent beyond about 4300 meters. 
Apparently this range defines the beginning of the asymptotic 
region for the conditions of these computer runs. Two things 
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are worth noting here: 1) the = const- similarity relation 

appears to work quite well, becoming better as the range 
increases, lending credence to our reasoning regarding the 
scaling o-f di -f -f r act i on beam spreading, whether arising from h or 
from R, and 2) remarkably, our result works well even tor 
propagation in the atmosphere, with its turbulence 
■f 1 uc tuat i ons (descr i bed as moderate for this case) and wind (but 
without thermal blooming). Hence, the similarity relation hold_. 
even beyond the range of validity implied bv the conditions of 
1 t s der 1 vat i on . 

Consul ting figure 4, in which the thermallly bloomer; i • e sno 
diameter is plotted, we see that adding the beam gu.a.lit / focroi 
to thie equal area circular aperture (curve 3) gets us only a 
little over half way to the 10"/- obscured aperture results 
2)- Hence, when the nonlinear effects of thermal biDornina are 
considered, the central obscuration is not so easilv compel . t -.3d 
for by our beam qu£-ility similarity reiatiori. However , tlio coce 
already handies the 107- obscuration case, ?.rid we are on ; \ 

[ ook {. nq a ineans of e;; tendi nq the r ai iq'-3 o . o;; i i i .a 

propagation codes to include larger obscurations and r ec t ancu i -c 
apertures. Clearl/, using the principle of constani 
part V4ay there- Hence, a possible approach is to qener t i a j 
e^xtend the beairi quality vs- R scaling relation '.or s i nn 1 _ r i t •/ 
i-elationi to the nonlirnear rGqion by a> : itting cui -es t*: 
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existing code output, and/or b) having identified the extent of 
the problem and a possible direction for a phenomenological fix, 
carry out studies with existing wave-optics codes on rectangular 
apertures and larger obscurations to obtain the applicable 
scaling laws. Possibly some gener al i zat i on of = const - 

which also accounts for beam power, e.g. R*M == f ^I) where 1 is 
the beam intensity, will prove adequate. 

III. AIM BIAS AND DISPERSION 

A number of the codes surveyed allow at most two options 
regarding focus of the laser beam. Either the user must assume 
E'Kact focus on the target, or the focus is at infinity. Wiieither 
this is justified or not should be determined for the specific, 
realistic pointing and tracking, and laser focusing systems 
being planned for. It is not difficult to show with some very 
simple cal cul at 1 ons that a range error, whether bias (constant 
from pulse to pulse.' or dispersion can lead to appreciable 
spreacJinq of the laser energy at the true target range- i-lo 
conclusion can be drawn without specific predictions for the 
rocusi nq system par arnet er s . 

With regard to bias of the laser system, a st r ai gh t f or war d 
analysis of the inputs and outputs from the HELhWS code indicate 
that there will be cases in which the low frequencv' atmospheric 
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turbulence and wind— induced aim bias, which, along with the 
postulated levels of laser system pointing bias will dominate 
the dispersion errors. For those cases with high levels at low 
frequency atmospheric turbulence and/or system aim bias, it mav 
be that the optimal firing will require either some type o+ 
pattern-f 1 r 1 ng of laser pulses or possibly an added dispersion 
(shot gun effect) to compensate. 

IV. CONCLUSIONS AND RECOMNENDAT IONS 

Based ofi our survey of a number of high energy laser 
pr opaqat 1 on codes, it would appear that none are niahiy p^■c_ls 
tools for evaluating Army laser weapons systems operatiiig . n 
realistic environments- That this mav be true is noc sun:n"i -.ir. 
for the wave optics codes which were generally desiqnea co be 
highly precise codes, valuable for baselining the iriore 
phenofnenol oq 1 cal scaling codes intended for systems level 
studies. However , they vjere usual! / not designed to include cl 
r a Pi go ot o+lects in an e-i f i c i eii t tiLunpul rilionai tor.w tc- L.*e u s 
1 1 1 s s 1 0 1 n s s t d 1 e s • 

However , consider tlie state/neni in L 5 J vot. I f.* , 1 . -o , f'.e-vj HK. 
beams may iiave r'ectanqular cr oss-sec 1 1 on s with nan -central rn 
non-c 1 r CLil ar obscurations with greater Kor less' than lo . 
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obscuration However, given the current primitive state 

of HEL weapon development and the intended systems anal yS 
applications of HELAWS, the BRLPRO assumptions for the beam 
Characteristics are sufficiently accurate." This may seriously 
underestimate the errors in approximating a rectangular aperture 
or one with larger than 107. obscuration (or both) with a 
suitably chosen circular aperture, possibly with 10/. 
obscuration. We feel that our discussion in section II oi this 
report is sufficient grounds to question the above quote. 

Hence, it may be advisable to consider some of the 
suggestions outlined in section II, above, reqaraing 
^modifications of present codes to more realistically handle 
rectangular apertures and larger obscurations. Ultimatelv f-he, e 
mr.ust be some validation either with a wave-optics code ^.such as 
the Livermore 4-D code, the NRL SSPARAMA code, etc) or with 
experiments on existing high energy laser systems. 
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